Fluctuation-dissipation relation and the Edwards entropy for a glassy granular 

compaction model 
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We analytically study a one dimensional compaction model in the glassy regime. Both correlation 
and response functions are calculated exactly in the evolving dense and low tapping strength limit, 
where the density relaxes in a 1/ In t fashion. The response and correlation functions turn out to 
be connected through a non-equilibrium generalisation of the fluctuation-dissipation theorem. The 
initial response in the average density to an increase in the tapping strength is shown to be negative, 
while on longer timescales it is shown to be positive. On short time scales the fluctuation-dissipation 
theorem governs the relation between correlation and response, and we show that such a relationship 
also exists for the slow degrees of freedom, albeit with a different temperature. The model is further 
studied within the statistical theory proposed by Edwards and co-workers, and the Edwards entropy 
is calculated in the large system limit. The fluctuations described by this approach turn out to match 
the fluctuations as calculated through the dynamical consideration. We believe this to be the first 
time these ideas have been analytically confirmed in a non-mean-field model. 
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Granular materials have had much experimental and 
theoretical attention in recent years. They are intriguing 
as they form an additional state of matter, fundamen- 
tally different from gases, liquids and solids 0. Specific 
for these systems is that the thermal energy of its con- 
stituents, i.e. the grains, is negligible compared to other 
relevant energy scales. Since the thermal energy is negli- 
gible, there is no inherent mechanism that makes a gran- 
ular system explore its phase space. Any energy fed to 
such system is quickly dissipated, and if left unperturbed 
the system becomes trapped in one of many metastable 
states with an essentially infinite lifetime. In many situ- 
ations where these materials are handled or used in pro- 
duction, they are continuously fed energy through ex- 
ternal perturbations. As a result the system starts to 
explore the available phase space and macroscopic quan- 
tities, such as the density, start to evolve. This situation 
has been experimentally examined |i2i,l3| through subject- 
ing a container filled with a granular material to many 
taps of acceleration F. The time between taps was large 
enough for the system to dissipate any excess energy, and 
settle into one of its many meta-stable states between 
each tap. The relaxation of density (or free- volume frac- 
tion) in these experiments was well fitted by an inverse 
logarithmic form. In a different context we have previ- 
ously 0, 1^ introduced a simple one dimensional model 
for which we were able to analytically derive this inverse 
logarithmic relaxation. We here consider a robust gener- 
alization of the model, and obtain both response and cor- 
relation functions through a direct dynamical approach. 
The system is further investigated in the context of Ed- 
wards' statistical theory of granular compaction 
and the result of the two approaches are compared. 

The model can be seen as a minimal model of the bot- 



tom layer in a granular material compacting under ta_ 
ping. It is related to the continuum car-parking model 
and consists of unit sized hard-core blocks (particles) 
positioned on a ring of length L. The blocks, which 
interact via hard-core repulsion, do a caged diffusion 
along the ring with diffusion constant D{V). The blocks 
are further able to evaporate from the ring at the rate 
rc(r) = exp(— /(r)). To model tapping induced diffusion 
and activated escape from the ring one could for example 
choose D{T) cx F and /(F) oc 1/F. When a gap of size 
larger than one opens up, we take it to be filled by a ran- 
dom deposition of a particle with the tapping-strength- 
independent rate — 0(F'^) > 0. This is meant to re- 
flect the fact that the gravitational pull on the particles 
is independent of the tapping strength. These rules are 
summarised in Figure ^ Since we are interested in the 
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FIG. 1: Schematic representation of the model with dynami- 
cal rules indicated. 



very slow dynamics exhibited by this model in the limit 
of a dense system subject to a weak tapping, we consider 
the weak tapping strength limit, in which we demand 
(consistent with diffusion and activated evaporation) 



re(F),i?(F)/e2 «rd, F \ 0, 



(1) 



where e is the free length fraction. It is this ordering of 
the rates, and not their actual form, that plays a crucial 
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role in determining the long time dynamics. The average 
time a gap of size larger than one, stays larger than one, 
before it is closed by diffusion, is proportional to 1/D{T). 
The time it takes to fill such a gap by a deposition event is 
proportional to 1/rd. Due to the above ordering of rates 
we see that in the weak tapping limit all gaps that open 
up will eventually be filled by a deposition. Therefore 
the effect on the long time dynamics of any evaporation 
events is suppressed. Thus we can use the effective rules 

To = 0, Td = oo, £'(r) — finite, 

for the long time evolution of the system. This shows 
that as long as the different rates satisfy Q]), then the 
long time dynamics is insensitive to the precise form of 
the evaporation and deposition rates. Taking this limit 
amounts to completely suppressing any fast processes, 
such as the evaporation of a block followed by a subse- 
quent deposition of a block in the created gap. 

In a different context we have previously jj, |5j derived 
the exact form of the density-density correlation function 
in the dense and low tapping strength limit for the spe- 
cific choice of D{T) oc T and /(F) oc l/F. This was done 
through a geometrical description of the problem, and we 
now extend the considerations to the present case, and 
include a calculation of the response function of the slow 
degrees of freedom. We view the time evolution of the 
system between deposition events as a diffusion of the 
gaps (between blocks) on the hyper surface of constant 
density 



N 



N 



A„ = i - iV, < A„ < 1 



(2) 



where Ajv is a vector consisting of all the N individ- 
ual gap sizes between N adjacent blocks. This is il- 
lustrated in Figure |21 for the case of a ring with only 
three blocks. Asymptotic (large system, low tapping 
strength, high density) dynamic properties of the sys- 
tem can now be calculated through assuming ergodicity 
on TTpf between deposition events, and using straightfor- 
ward geometrical considerations. Ergodicity between de- 
position events should hold in the dense limit since the 
time between deposition events diverges, and thus the 
diffusive motion has enough time to relax the system be- 
tween these events. For further technical details we refer 
to nil [13 ■ (The intermediate time regime, in which the 
system ought to display spatial structures, is at present 
being investigated [ll|). The connected two time density- 
density correlation function has the form 



T2{ew, £o) et 
LT(e)T(ew) ^ 2ie2 



=.l/cw-l/e 



with e = e(t|eo,io) and ew ~ e(tw|eo,io), where e(i|eo,io) 
denotes the average free-length fraction on the ring at 
time t, given the initial free-length fraction eg at time 
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FIG. 2; Difi^usion on the hyperplane showing a few reflections 
at the boundary, corresponding to two blocks bouncing off 
each other, and an eventual escape and transfer to the hyper 
plane in through a gap of size 1 opening up and then being 
flUed with a block. 



Iq. In the above is the waiting time, and T(e) ~ 
ke^e^^^ / DiV) is the average time between deposition 
events within a unit length of the ring, given the free- 
length fraction e. Here fc is a constant only dependent on 
the geometry of the hyper planes. We have further used 

T„(ew,eo) = y" °der"(e) 
k 



D{T) 



i(e^("+i)(l + 0(ew))e"/^" 
-ef'+^^(l + 0(6o))e"/^° 



Thus we see that the long time density-density correla- 
tion is independent not only of the evaporation and de- 
position rates, but also of the diffusion constant. Since 
Ti(e,eo) by definition is the average time it takes the 
system to evolve from a free-length fraction eo to a free- 
length fraction e, we have an implicit relationship for the 
evolution of the free-length fraction 

Ti(e(t|eo,to),eo) ^ t - to. 

Through differentiating the above with respect to F it is 
an easy matter to calculate the response in density to a 
change of the tapping strength starting at time tv,: 



Xr(t,iw) 



c>e(^|ew,^w) _ D'{r)Ti{e,e^) 
dT 
e^D'jT) 
D{r) 



1 



D{r)r{e) 
(^)%xp(l/e„-l/e)^ 



In Figureiniwe display the response as a function of time. 
Since we have both response and correlation functions 



3 



Xr(t,tw) 



10** 



10^ 



FIG. 3: The form (scale arbitrary) of the response function of 
the slow degrees of freedom, xr(i, iw), as a function of t — tw, 
and with a initial free-length fraction at tw, tw = 0.05. Similar 
response functions have previously been found in 0, (l^J . 



we are in position to consider a possible extension of the 
equilibrium fluctuation-dissipation theorem to the non- 
equilibrium situation present in our model. A direct 
comparison of the asymptotic form for the response and 
connected correlation functions gives 

2 T D'(V) 

Xr(t, iw) {C^{t, t) - C7^(i, t^)) . 

In an equilibrium system the fluctuation-dissipation the- 
orem states that the linear response and the correlation 
function (for intensive quantities) are related through 



cq 



where V is the system size and h is the variable conjugate 
to the quantity considered. In our case the corresponding 
statement is 



Tr 



ncq 



h'{r)D{T) 

2D'{r) 



where h(T) is the (unknown) variable conjugate to the 
density. It should be noted that in our asymptotic analy- 
sis we have nowhere demanded that i S> iw, but just that 
t — t^ is larger than the relaxation time of the fast degrees 
of freedom, ttast- Thus the above results are valid as long 
as S> to and t—t^^ ifast- Therefore Tneq could in gen- 
eral depend on some finite combination of t and iw in the 
limit — > oo- It is thus striking that Tnoq is completely 
independent of t and t^. The temperature of the fast de- 
grees of freedom can be directly calculated ilfl] through 
the fluctuation-dissipation theorem as Teq = h' / f , and 



therefore Tcq/Tn, 



2D'/{Df'). Thus, depending on 



which timescales we are considering, we see different dy- 
namically defined temperatures. This very behavior has 
previously been identified in mean-field models [l^ and 
through numerical simulations 0, Hil IT^ , and now we 
see it analytically in a non-mean-field model. For any rea- 
sonable system we have D' > and /' < 0, and thus the 
two temperatures differ in sign (which has not been seen 
in mean field models). This sign difference arises because 
on the short timescale an increase in tapping strength 
decompactifies the system through the fast degrees of 



freedom considered above. In the aging regime though, 
a higher tapping strength increases the compactification 
rate. This behavior has previously been seen in systems 
with activated dynamics, such as the trap model (see 
and references therein). We will later compare these re- 
sults with those from considering the Edwards entropy 
of the system, which we derive next. 

We will use the simplest form of this theory and apply 
it to our model, and for the omitted technical details we 
refer to 0,^3- ^ similar treatment has recently been car- 
ried out for the parking-lot model jl9| , and similar results 
for the simpler case when one allows gaps of size larger 
then one has been known for a long time [20j . The count- 
ing of the number of blocked configurations is central to 
the Edwards approach, and in our case these correspond 
to configurations for which no gap is larger than one. A 
powder driven by well separated periodic taps will ex- 
plore the phase space of meta-stable or blocked states, 
consistent with the external conditions. It is then natu- 
ral to define the entropy density in the thermodynamic 
limit as 

SEdw(e) = lim ^\nWL{e), 

where Wl (e) is the number of blocked microscopic states 
consistent with the free-length fraction e. We will refer 
to SEdw as the Edwards entropy density. In analogy with 
equilibrium statistical mechanics one assumes that any of 
the states consistent with the macroscopic constraints e 
and L are equally probable. With this crucial assumption 
the statistical properties of the system are given by the 
micro-canonical partition function 

Z(e, L) = exp(LsEdw(e))- 

If we consider our original system as being part of a larger 
ensemble that allows exchange of particles between its 
subsystems, then we can move over to the canonical en- 
semble. We define the canonical partition function for 
the free-length fraction by 



n{^,L)~ / deexp(— L(/ie 

^0 



— SEdw(e)), 



where /i is a Lagrangian multiplier ensuring the correct 
overall free length. In the thermodynamic limit we can 
use the saddle-point method to introduce the thermody- 
namic potential 

/(m) = - lim y \nn{^, L) = fie{fi) - SEdw(e(M)) 

where e(/i) solves 5eSEdw(e(M)) = In the usual man- 
ner we can now calculate the average free-length fraction 
as a function of /j,, (e)^ = 9/i/(A«) = e(A<). Through the 
assumption that all possible states have a priori equal 
probability we have thus, with respect to the calculation 
of the steady state, been able to cut out any reference to 
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the dynamics and replaced it with a statistical descrip- 
tion using the thermodynamic variable /x conjugate to 
the free-length fraction. In our present setting the a pri- 
ori equal probability is just the earlier used assumption 
of ergodicity on ttjv, which should hold true in the dense 
limit. The Edwards entropy for the compaction model 
considered can then simply be expressed as 

SEdw(e) = lim ^ lnVol(7r2,(i_e)). 

L—i-oo L/ 

where Vol(7r2,(i_e)) is the volume of the hyper surface 
7ri(i_e) defined in (jSJ. Asymptotic analysis j^, E| gives 

SEdw(e) = -g(e,s(e)) 
g{e,s) = 2es — (1 — e)(lnsinhs — Ins + s), 

where s(e) solves dsg{e, s{e)) = 0. Moving over to the 
canonical ensemble we need to solve /i = 9eSEdw(e(M)) for 
the average free-length fraction (e)^ = e(/i) as a function 
of the thermodynamic variable /i. In the low tapping 
strength, high density, and large system limit this can be 
calculated as 

SEd„(e)-^(l-e)(l-ln(e-i-l)) 

and thus e(/i) = + 0(m^^ 1^^/^))- The density fluc- 

tuations can be written as 

(oe )Edw = --^Om^IW ^ 

which should be compared with the fluctuations of the 
slow degrees of freedom, C'^{t,t) ^ e^/2L, as given by 
the dynamical considerations above. Though different 
in origin they agree in their functional dependence on e. 
Results with similar implications have previously been 
reported in numerical simulations of compaction mod- 
els [13 El El. 

Lastly we note that in a steady state the free-length 
fraction, e(r), is set by the condition that the effective 
evaporation rate for a typical configuration matches the 
deposition rate into gaps opened by diffusion, 

exp(-(/(r) + r^e\r)/D{T))) ^ exp(-l/6(r)). 

For the natural choice f{T) = 1/kT, and D{T) oc T, this 
gives e(r) ~ kT, which compared to e(/x), as given by the 
Edwards picture, gives /i 1/kT. Thus, in this setting 
it is natural to interpret F as a temperature for the slow 
degrees of freedom. 

In this letter we have presented what we believe is 
the first instance where an asymptotically exact calcu- 
lation gives an extension of the fluctuation-dissipation 
theorem in the non-equilibrium regime of a granular com- 
paction model. The results are very robust with respect 
to the different types of driving, with only the relation 
setting the limits for the possible forms. Through the 
fluctuation-dissipation theorem it is possible to define 



two different temperatures, one for the fast and one for 
the slow degrees of freedom. Due to the difference in the 
response to a perturbation in the tapping strength of the 
short time and long time degrees of freedom, these tem- 
peratures have a different sign. We have further treated 
the model with the Edwards theory of powders and cal- 
culated an exact (implicit) form of the Edwards entropy 
density. The fluctuations as calculated within the Ed- 
wards' picture further accurately describe the long-time 
dynamically induced fluctuations. 

The authors gratefully acknowledge Juan P. Garrahan 
and an anonymous referee for valuable comments. This 
work was supported by EPSRC under the Oxford Con- 
densed Matter Theory Grant No. GR/R83712/01 and 
GR/M04426. MD gratefully acknowledges support from 
the Merton College Oxford Domus scholarship fund, and 
the Royal Swedish Academy of Science. 



* MD's present address: Instituut-Lorentz, Leiden Univer- 
sity, P.O. Box 9506, 2300 RA Leiden, Netherlands 
email: depken@lorentz.leidenuniv.nl 

[1] H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Rev. 
Mod. Phys. 68, 1259 (1996). 

[2] J. B. Knight, C. G. Fandrich, C. N. Lau, H. M. Jaeger, 
and S. R. Nagel, Phys. Rev. E 51, 3957 (1995). 

[3] E. R. Nowak, J. B. Knight, E. Ben-Nairn, H. M. Jaeger, 
and S. R. Nagel, Phys. Rev. E 57, 1971 (1998). 

[4] R. Stinchcombe and M. Depken, Phys. Rev. Lett. 88, 
125701 (2002). 

[5] M. Depken, Ph.D. thesis. Department of Physics, Uni- 
versity of Oxford (2003). 

[6] S. F. Edwards and A. Mehta, J. Phys. France 50, 2489 
(1989). 

[7] A. Mehta and S. F. Edwards, Physica A 157, 1091 
(1989). 

[8] R. B. S. Oakshot and S. F. Edwards, Physica A 189, 188 
(1992). 

[9] A. J. Kolan, E. R. Nowak, and A. V. Tkachenko, Phys. 

Rev. E 59, 3094 (1999). 
[10] M. Depken and R. Stinchcombe (2004), to be published. 
[11] A. Lefevre, L. Berthier, and R. Stinchcombe (2004), to 

be bublished. 

[12] A. Barrat and V. Loreto, J. Phys. A 33, 4401 (2000). 
[13] M. Nicodemi and A. Coniglio, Phys. Rev. Lett. 82, 916 
(1999). 

[14] L. F. Cugliandolo and J. Kurchan, J. Phys. A 27, 5749 
(1994). 

[15] G. Parisi, Phys. Rev. Lett. 79, 3660 (1997). 

[16] J.-L. Barrat and W. Kob, Europhys. Lett. 46, 637 (1999). 

[17] A. Barrat, J. Kurchan, V. Loreto, and M. Sellitto, Phys. 

Rev. Lett. 85, 5034 (2000). 
[18] P. SoUich, S. Fielding, and P. Mayer, J. Phys. : Condens. 

Matter 14, 1683 (2002). 
[19] G. Tarjus and P. Viot, Phys. Rev. E 69, 011307 (2004). 
[20] L. Tonks, Phys. Rev. 50, 955 (1936). 
[21] A. Barrat, J. Kurchan, V. Loreto, and M. Sellitto, Phys. 

Rev. E 63, 051301 (2001). 
[22] M. Sellitto, Eur. Phys. J. B 4, 135 (1998). 



